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Abstract. In the area of traditional physics the atomic nucleus belongs to the most complex sys- 
tems. It involves essentially all elements that characterize complexity including the most distinctive 
one whose essence is a permanent coexistence of coherent patterns and of randomness. From a more 
interdisciplinary perspective, these are the financial markets that represent an extreme complexity. 
Here, based on the matrix formalism, we set some parallels between several characteristics of com- 
plexity in the above two systems. We, in particular, refer to the concept - historically originating 
from nuclear physics considerations - of the random matrix theory and demonstrate its utility in 
quantifying characteristics of the coexistence of chaos and collectivity also for the financial mar- 
kets. In this later case we show examples that illustrate mapping of the matrix formulation into the 
concepts originating from the graph theory. Finally, attention is drawn to some novel aspects of the 
financial coherence which opens room for speculation if analogous effects can be detected in the 
atomic nuclei or in other strongly interacting Fermi systems. 
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COMPLEXITY 

The concept of complexity refers to diversity of forms, to emergence of coherent patterns 
out of randomness and also to an extreme easiness of switching from one such pattern 
to another. Such properties typically accompany the systems that involve many degrees 
of freedom, many different space and time scales, and therefore such phenomena like 
chaos, noise, but also coherence, collectivity and criticality [1]. It thus seems most ap- 
propriate to view the real complexity just at the interface of chaos and collectivity [2, 3] 
as these two seemingly contradictory phenomena involve essentially the same elements. 
The recent concept of the scale free complex networks [4] offers one promising direction 
to formalize complexity though many related issues still remain open. 

Studying complex systems, either empirically or theoretically, is unavoidably based 
on analyzing large multivariate ensembles of parameters. The most efficient starting for- 
malism to grasp such ensembles and to quantify the whole variety of effects connected 
with complexity is in terms of matrices [5]. Since complexity is embedded in chaos, or 
even noise, the random matrix theory (RMT) [6, 7] provides then an appropriate ref- 
erence. Its utility, at the first place, results from the fact that the degree of agreement 
quantifies the generic properties of a system - those connected with an overwhelming 
chaotic or noisy dynamics. For the complex systems this is expected to be a dominant 
component, but this component is not what constitutes an essence of complexity. From 



this perspective the deviations are even more relevant and more interesting as they may 
reflect constructive and perhaps deterministic components emerging from a noisy back- 
ground. 

The main related purpose of the present summary is to explore the financial char- 
acteristics of complexity in analogy to the ones that already are identified in strongly 
interacting quantum many body systems such as atomic nuclei. One such principal char- 
acteristics in the nuclear context is the energy separation of the collective state from the 
noisy background [5, 8, 9]. 

COMPLEX NUCLEI 

The structure of nuclei is a conventional example for collectivity embedded in chaos 
[8, 9]. The single particle structure in odd mass nuclei, the rotational spectra in deformed 
nuclei, the low-lying collective 2+ and 3~ states in spherical nuclei as well as the high- 
lying giant resonances are indications of a regular behavior in nuclei. The absolute 
number of such states in each nucleus, however, is small compared with the total number 
of levels. For that reason statistical methods have been developed and used to analyze 
the remaining huge number of individual levels [10]. The best example for collectivity 
embedded in chaos are giant resonances. These collective states exist in all medium 
and heavy mass nuclei. In Figure 1 the electric giant quadrupole resonance in several 
nuclei is shown. Collectivity means a cooperation, and thus the coupling between the 
different degrees of freedom in order to generate a coherent signal in response to an 
external perturbation. In shell model approaches collective states may be generated 
within the random phase approximation (RPA) in a lplh configuration space [11]. For 
a quantitative description of the width of giant resonances, however, it is necessary 
to couple to the 2p2h space [12] where, due to a sufficiently large density of states 
relative to the strength of the residual interaction, local level fluctuations characteristic 
for a Gaussian orthogonal ensemble (GOE) appear to already take place [8, 13]. The 
collective lplh part of the giant resonances can thus be considered as embedded in the 
chaotic 2p2h space. The mechanism that creates coherent solutions can transparently 
be understood in a schematic model postulated by Brown and Bolsterli [14]. Here 
one assumes that all lplh states are energetically nearly degenerate as indicated in the 
upper part of the right panel of Figure 1 . As interaction one uses a separable multipole- 
multipole interaction which allows one to solve the RPA equation analytically. The result 
is shown in the lower part of the right panel of Figure 1 . One state is strongly shifted 
in energy and includes all the transition strength B(Q). All the other states loose nearly 
all their transition strength. In a more realistic description, if the 2p2h space is properly 
included, one is able to explain the experimental spectra shown in Figure 1 [12]. 
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FIGURE 1. Example from the nuclear physics. (Left) Giant quadrupole resonances in atomic nuclei. 
(Right) RPA solution in the schematic Brown-Bolsterli model; in the upper part the degenerate lplh states 
are shown, in the lower part the corresponding RPA solution for a repulsive ph-interaction. 



COHERENCE VERSUS RANDOMNESS IN THE MATRIX 

REPRESENTATION 



Expressed in the most general form, in essentially all the cases of practical interest, the 
nxn matrices W used to describe the complex system can be written as 



W = XY T , 



(1) 



where X and Y denote the rectangular n x m matrices. Such, for instance, are the 
correlation matrices whose standard form corresponds to Y = X. In this case one deals 
with n observations or cases, each represented by an m dimensional row vector x ( - 
(y,), (z = l,...,n) [15], and typically m is larger than n. In the limit of purely random 
correlations the matrix W is then said to be a Wishart matrix [16]. The resulting density 
pw(A) of eigenvalues is here known analytically [17, 18], with the limits (A, ran < A < 



X max ) prescribed by V™* = 1 + 1/(2 ± 2 y^l/Q and Q = m/n > 1 : 
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This is the so-called Marcenko-Pastur distribution. The variance o of the elements of x,- 
is in our case equal to unity. 

As an appropriate reference for quantifying various characteristics of complexity one 
can use a few distinct types of random matrix ensembles, depending on a problem. In 
a typical situation, the bulk of such characteristics is consistent with RMT predictions. 
However, the genuine properties of a studied system are expected to manifest themselves 



in deviations from the universal pattern. These deviations are primarily connected with 
the reduction of dimensionality due to emergence of synchronous or coherent patterns. 
In the matrix formalism such patterns can be expressed with a significantly reduced rank 
of a leading component (or components) of W. This can be written schematically as a 
decomposition: 

W = W RMT + W coll + £W i -, (3) 

in which the first component is the noisy (RMT) component, the second (usually dom- 
inating) component of rank 1 describes the coexisting collective global behavior of the 
system, while some possible additional components denoted by the sum are related to 
local coherent effects engaging only parts of the system. 

In the absence of Wrmt, the other terms in (3) generate k nonzero eigenvalues, and 
all the remaining ones (n — k) constitute the zero modes. The schematic Brown-Bolsterli 
model case discussed above corresponds to k = 1 with no other terms in (3). When - more 
realistically - Wrmt enters as a noise correction, a trace of the above effect is expected 
to remain, i.e., k large eigenvalues and the bulk composed of n — k small eigenvalues 
whose distribution and fluctuations are (of course after disregarding the collective part) 
consistent with an appropriate random matrix ensemble. A likely mechanism that may 
lead to such eigenspectra is that m in Eq. (1) is actually significantly smaller than n, or 
that it is effectively small if the number of large X and Y entries is small. 

Another mechanism that may lead to a structure analogous to (3) is the presence 
- in addition to noise - of some systematic trend in the column vectors of X and Y. 
Then [5, 19], to a first approximation, the W co n term is represented just by a matrix 
whose all entries are close in magnitude, and thus its rank is directly seen to be unity. 
This structure of a matrix can be easily detected if one observes an asymmetric shift of 
P(w) relative to zero. 

The matrix form as expressed by Eqs. (1) and (3) proves very general. The Hamilto- 
nian matrices of strongly interacting quantum many body systems such as atomic nuclei 
can also be expressed in this form. This holds true [20] on the level of bound states where 
the problem is described by the Hermitian matrices, as well as for excitations embedded 
in the continuum [21]. 

FINANCIAL MARKETS 

Financial markets have only relatively recently become area of non-episodic scientific 
interest for physicists. Despite some interesting earlier connections between physics 
and financial economics like the description of the price movements in terms of the 
random walks [22] or the approach to modeling of bond prices by means of the path 
integrals [23], physics and economics for decades were developing largely indepen- 
dently. This changed rapidly after financial institutions absorbed a number of physicists 
in 1990's and, at the same time, large amounts of computer-recorded financial data be- 
came easily available. Financial data are alluringly clean and this forms an additional 
factor attracting physicists and, among them, especially the statistical physicists. 

From an empirical perspective, financial markets, like the stock markets or the cur- 
rency markets, are systems with many degrees of freedom defined by the assets traded 




FIGURE 2. Example from the financial physics. (Left) Distribution of correlation matrix elements for 
time series of 15 min. returns of 1000 American companies traded on NYSE and Nasdaq. Inset: the same 
but in the double linear scale. (Right) Eigenvalue spacings distribution for the same correlation matrix 
(histogram) compared with the Wigner distribution predicted for GOE matrices by RMT (dashed line). 



on a given market. Since price movements are inevitably dominated by trends that make 
data strongly nonstationary, it has become standard to analyze time series of price re- 
turns (defined as logarithmic price increments over fixed time intervals) instead of prices 
themselves. It does not solve the problem of nonstationarity completely, but at least al- 
lows one to successfully use the standard methods of time series analysis. 

In accord with what has been presented in the previous sections, also a financial 
market can be described as a system of interacting or otherwise coupled elements-assets 
whose dependences can be expressed by a matrix. In the present case this role is played 
by a correlation matrix which, in accordance with Eq. (1), is constructed from a data 
matrix X consisting of multivariate time series of returns of n assets observed over 
equivalent time intervals. After diagonalization of the correlation matrix, one obtains 
its eigenspectrum which then can be compared with the RMT prediction for the Wishart 
ensemble. 

As an illustration, we consider a set of time series of 15-min. returns for n = 1000 
highly-capitalized American companies traded on the New York Stock Exchange or on 
Nasdaq. Each time series is m = 13,546 points long and comprises the interval 1998- 
1999. First of all, the left panel of Figure 2 shows the distribution of matrix elements 
which is the most straightforward measure of the matrix properties. The distribution is 
asymmetric and, exactly as discussed above, its center of mass is shifted towards positive 
values. As this Figure documents, the distribution is also far from Gaussianity. 

What is now completely obvious, the eigenspectrum of our correlation matrix also 
deviates significantly from the RMT result [24]. This maybe is not so well-evident on the 
level of the eigenvalue spacings (right panel of Figure 2) whose distribution is roughly 
consistent with the Wigner distribution (the Wishart random ensembles do not differ in 
this respect from the Gaussian orthogonal ones). However, if one directly looks at the 
eigenvalue spectrum (left panel of Figure 3), there is very clear non-random structure 




FIGURE 3. (Left) Eigenvalues A,- of the correlation matrix calculated for time series of 15 min. returns 
of 1000 American companies traded on NYSE and Nasdaq. Note the logarithmic scale of the axis. 
(Right) Eigenvalue density for the same matrix (histogram) together with the corresponding Marcenko- 
Pastur distribution predicted by RMT (dashed line). This empirical distribution does not show the four 
largest eigenvalues which, as it is seen in the left panel, are located beyond the range of the graph. 
Inset: Marcenko-Pastur distribution in agreement with the eigenvalue distribution for a correlation matrix 
constructed from the randomized data (histogram). 



of especially the largest several A's. The largest eigenvalue Ai = 84.3 is extremely 
strongly repelled from the rest of the spectrum. This eigenvalue can be identified with 
the W C oii term in Eq. (3), which in the language of the financial markets is called "the 
market factor" and describes the collective movement of all assets (usually grasped by 
the behavior of a market composite index). A few smaller eigenvalues seen in Figure 3 
have to be associated with the third term in Eq. (3); they typically correspond to market 
sectors (industries in the case of stock markets or geographical regions in the case of the 
currency market). It must be noted that the RMT bounds are quite narrow in the present 
case: A^ = 0.53 and Amax = 1.62. 

A characteristics which allows one to quickly assess the existence of non-random 
components in the correlation matrix is the eigenvalue density. This is because of 
a concise result for the Wishart matrix ensemble expressed by (2). This theoretical 
distribution is shown in the right panel of Figure 3 together with the empirical eigenvalue 
distribution for the data under study. In fact, the discrepancies between the two are 
significant both at the ends and in the center. The eigenvalue density for the randomized 
time series (surrogates) is seen in the inset. This proves that the non-random properties 
of the actual eigenvalue density come from the inter- stock correlations destroyed by the 
shuffling of data. 

The results presented so far describe the statistical properties of dependences between 



FIGURE 4. (Left) Minimal spanning tree for a set of 1000 large American companies traded in NYSE 
and Nasdaq. (Right) MST for a subset of 100 American companies of the largest capitalization; edge 
widths are proportional to the corresponding metric distances (4). 

the assets but do not allow us to visualize the exact microscopic structure of these 
dependences. This may be considered a drawback especially in the case of financial 
markets where, for practical reasons, the details are important. This issue can, however, 
be overcome with the network-oriented approach, in which the relevant relations among 
the system's elements can be expressed by a structure of nodes and edges. In order to 
illustrate the point, we create a network representation of our set of 1000 American 
stocks. Each asset is defined as a node and each pair of assets is connected by a link with 
a weight equal to the corresponding correlation matrix entry. However, the so-defined 
network consists of 1000 nodes and 495,000 edges which makes it virtually impossible 
to be presented graphically. Therefore we have to confine our task to the construction of 
a sub-network of "minimum weight", i.e. the sub-network in which connections between 
nodes of each pair are the strongest possible. This is the so-called minimum spanning 
tree (MST). 

The MST is calculated based on the correlation matrix after transforming its entries 
Wij into metric distances [25] according to the following relation: 



The distances dn are sorted from the smallest to the largest forming a list and by going 
from the top of this list one connects pairs of nodes which have not yet been connected 
to each other by any continuous path. In this way, after passing through the whole list, 
one finally obtains a dendrogram with n nodes and n — 1 edges, which is relatively easy 
to draw and analyze. MST for our set of assets is displayed in Figure 4 (left panel). It is 
striking that the resulting network consists of a central group of few stocks to which all 
other stocks, in a form of layers, are connected. This can be interpreted as an existence of 




(4) 



companies which are treated by investors as the whole American economy in a nutshell. 
Other companies are thus viewed as peripheral and less representative. 

In order to better see the core structure of the market and identify the "hubs", the 
same procedure was applied to the subset of 100 companies with the largest market 
capitalization (roughly > 3 • 10 10 USD). The corresponding MST is shown in the right 
panel of Figure 4. Now it is clear that the most important node is General Electric (GE) 
and it is accompanied by less central, but still important, nodes of Citigroup (C) and 
Cisco Systems (CSCO). GE is a huge conglomerate with diversified areas of activity and 
this fully justifies its dominant role in the network. On the other hand, C is a part of a 
large industry of financial institutions and CSCO belongs to the high-tech market sector 
which is also significant. Other companies which share the entire or a part of activity 
areas of the dominant companies (or are otherwise related) are naturally connected to 
them by direct links, forming the inner layers. In the same manner, the rest of nodes 
form the outer layers since their profiles significantly differ from the profiles of the three 
hubs. 

This example of MSTs shows that if one needs to grip the key features of the micro- 
scopic structure of a system quickly and in a concise manner, the network formalism 
proves very useful [26, 27]. But also more general and global characteristics of de- 
pendences among the system's constituents can be described using approaches based 
on the networks theory. This refers primarily to the case of complex systems in which 
complexity manifests itself by such features like scale invariance or the small world 
structure [28, 29] and the matrix representation offers an appropriate starting point. 

SUMMARY 

As shown above, several concepts originally developed in the area of nuclear physics and 
the related experience of handling complexity prove useful in attacking even more com- 
plex systems such as the financial markets. Some global characteristics like the mech- 
anism of the coexistence of chaos and collectivity appear to be similar. As far as the 
financial collectivity is concerned there are some novel elements however. One of the 
related distinct patterns, likely originating from the hierarchical organization illustrated 
in the previous Section, are imprints of the discrete scale invariance which manifest 
themselves via the log-periodic oscillations [30] cascading self-similarly through vari- 
ous time scales [31] with amazingly universal quantitative characteristics [32, 33]. It is 
not excluded that at a sufficiently high energy resolution the analogous effects can be 
detected also in the atomic nuclei. Finally, the insight offered by the network represen- 
tation, as above illustrated for the stock market, may encourage attempts to graphically 
represent the nuclear Hamiltonian as well, in the mean field basis for instance. 
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